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THE PERMEABILITY OCF POROUS SOLIDS TO 
GASES AND LIQUIDS.* 


By L. Grunsere, B.Sc., and A. H. Nissan, Ph.D. 


* INTRODUCTION. 


CALCULATIONS reveal that the rates at which fluids flow through sand 
formations are extremely small. Thus an oil-well flowing, say, at 2000 
tons/day from a sand 30 feet thick may be considered to be a fairly prolific 

ucer. Yet, if the specific gravity of the oil is 0-8, the mean linear 
velocity at a distance of only 3 feet from the centre of the well is of the 
order of 0-05 cm./sec. The speed of flow is inversely proportional to distance 
from the well and, hence, it can be seen that in wells of more normal size 
the mean linear speed of the fluid in the body of the sand quickly assumes 
exceedingly small values. In laboratory work, however, there are 
advantages in maintaining fairly high rates of flow as the measurements 
become relatively simpler to take and the accuracy easier to maintain. 
With these considerations in mind, it was deemed necessary to investigate 
the validity of certain principles connected with the flow of fluids through 
porous media, when the rates of flow become very low. “To attain easily 
measured pressure drdp across the samples at low rates of flow of air 
through them, highly impervious samples were chosen. Liquids were 
excluded from the studies on rock samples for two reasons: (1) liquids 
introduce complex extraneous factors, such a3 the swelling of certain 
minerals, which factors may mask fundamentai deviations in the hydro- 
dynamics of the system; (2) certain measurements have already been 
taken by earlier workers. To obtain some relationship between the flow 
of air and the flow of liquids a Jena glass filter was used. In this case the 
difficulty mentioned above was eliminated. 

Permeability is the property of porous materials, which permits the flow 
of liquids through these materials. A porous material is permeable by 
virtue of its porosity, but any attempt to express permeability as a function 
of porosity has failed so far. Clays, for example, having a porosity of 
50-60 per cent. are practically impermeable, while sandstones with a 
porosity of 30 per cent. are quite permeable. However, it can be safely 
stated that permeability is a function of the dimensions and intercom- . 
munication of the pores, factors which in most cases cannot be determined. 


Theory. 

The fundamental principle governing the flow of fluids through porous 
media, when the flow is in the viscous region, is “ Darcy’s law.”+ This 
law in its simplest form states : 


dp 


* Received October 30th, 1941. 
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where yd = pressure gradient along sample ; 
quantity passing/unit 
cross-sectional area 


v = macroscopic velocity = 
C = a constant. 


The constant C contains a physical characteristic of the fluid, i.c., its vis. 
cosity u. It ig usually assumed that separation of u from the constant ¢ 
frees the latter from any dependence on the nature of the fluid. In a more 
elaborate form Darcy’s law becomes 
Qu ii 

where K = permeability—i.c., a characteristic of the porous medium; 

Q = volumetric of flow with respect to time ; 

= viscosity of the fluid ; 

A = cross-sectional area of the sample ; 

AP = pressure drop per unit length. 


‘In the case of flow of gases the quantity Q is not constant, but increases 
with the pressure drop, according to Boyle’s law, in the majority of cases, 
It is therefore necessary to correct this quantity, as measured at the outlet, 
to the mean pressure. Assuming isothermal expansion : 


(3) 
where Q = mean quantity flowing per unit time ; 
P, + P, = absolute pressures at inlet and outlet respectively. 


For gases, therefore, Darcy’s law becomes : 


@ 


It must, however, be pointed out that Darcy’s law is, in fact, another 
expression of Poisseuille’s law of flow through a number of capillaries : 


(5) 


where n = number of capillaries ; 
R = average radius of capillaries. 


nR® 


Therefore K=C’ 


where C” is a constant introduced by the fact that in Darcy’s law, not the 
actual cross-sectional area of the capillaries through which flow is taking 
place, is used, but the macroscopic cross-sectional area (i.e., the cross- 
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sectional area of the sample) and further, the cross-section of the capillaries 
js neither uniform nor circular. 

This relationship between Darcy’s and Poisseuille’s laws will later be 
ysed to find the correct values of certain criteria which were found to be 
controlling factors of the conditions of flow. 

In the c.g.s. system K is measured in darcies, when Q (or Q) is measured 
ine.c./sec., A in cm.?*, » in centipoises, and AP in atmosphere per centimetre 
length. The practical unit is the millidarey = 10 darcy. 

The absolute validity of Darcy’s law was challenged by early workers, 


who found that the quantity v did not vary linearily with =. but varied 


according to an exponential law. One form * proposed was : 


(7) 


where a and n are constants. 
Another form given to this equation was : 


where a, 6, and ¢ are constants. 
Based on an analogy with the flow of fluids through pipes, the following 
form is far the most popular :* ? 
dp 

dx 


=av+b®. . ..... (9) 


where a and 6 are constants. 

Based on the above analogy and on mathematical analysis, it is believed 
that for low velocities the fluid is in a state of “ viscous flow,” governed 
almost entirely by the first term of equation 9. As the velocity increases the 
fluid changes to a state of “‘ turbulence ’’ governed by the second term of 
equation 9. M. Muskat limits the importance of the linear term to systems 
with Reynolds’ criterion lower than 1. The parameter for length in the 
criterion has hitherto been difficult to determine, so that this figure may not 
represent any strict limitation of the conditions under which viscous flow 
changes to turbulence. 


(4) 
other 


(5) 


Oxssect OF INVESTIGATION AND CHOICE OF MATERIAL. 


As has already been mentioned, it was the object of the following investi- 
gation to test the validity of Darcy’s law for rock samples, which, due to 
their small pore “diameter, will, under test conditions, give non- turbulent 
flow. 

Oolitic limestone was chosen, which is known to have approximately 
73 per cent. of its total porosity in the form of micropores, with an equivalent 
diameter smaller than 0-005 mm. The results were to be correlated by tests 
on limestone, magnesium limestone, fine (siltstone-like) sandstone, and 
coarse (well-cemented) sandstone. The impervious samples were used to 
yield extremely low rates of flow, whilst the relatively permeable sandstones 
were used to correlate low and high rates in one system. 
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A Jena glass filter (commercial specification 11G3) was used to correlate 
the flow of air and of aqueous solutions and to study the permeability to 
liquids of different surface tensions. 


EXPERIMENTAL WoRK ON Gas FLow. 


The cores were cut by means of a special cutter, and whenever necessary 
subjected to Soxhlet extraction, according to a standard procedure.* The 
dimensions of each sample were taken by means of a micrometer graduated 
in millimetres. Length and diameter were measured in five different 
planes, and the mean of five readings taken. The results are given in 
Table I. 


Taste I. 


Lithological species. | Length, cm. | Diameter, cm. aon 


Oolite 3-1936 2-9533 6-8500 
Oolite 3-1991 2-9656 6-9074 
Oolite 3-1916 2-9528 6-8478 
Oolite 3-1996 2-9728 6-9411 
Limestone | 3-2536 2-8589 6-4192 
Very fine sandstone | 1-9684 3-1003 7-5491 
Mg. limestone 3-0304 2-8390 6-3258 
Coarse sandstone 2-3879 2-9563 6-8799 


Coarse sandstone 2-3499 2-9668 6-9130 


The samples were thoroughly dried for 6 hours at about 70° C., and 
always kept in a desiccator over calcium chloride. The sample was placed 
in the special sample holder, and air under pressure passed through. Keep.- 
ing the pressure constant, measurements were taken simultaneously for: 
barometric pressure, temperature, inlet and outlet pressures and rate of 
flow. 

Special precautions were taken to keep the sample in a state of dryness 
throughout the experiment, firstly by passing the air through a calcium 
chloride container, and secondly by interposition of calcium chloride tubes 
whenever there was communication to measuring instruments containing 
water. 


APPARATUS. 
The apparatus (Fig. 1) may be sub-divided into three sections : 


(1) Source of pressure and sample holder ; 
(2) Pressure-measuring section ; 
(3) Flow-measuring section. 


(1) Air under pressure from a compressor passes through two reducer 
valves (RV, and RV,) and two pressure tanks (P7', and PT-,) into a calcium 
chloride container (CC,). The latter also holds the thermometer (7h), 
which measures the temperature of the stream of air. Passing through a 
filter of cotton-wool, the air goes through a T-piece to the sample holder 
and to the pressure-measuring section. 
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The sample holder (SH and Fig. 14) consists of a steel container and 
different brass tubes and rings. Rubber packings are provided to prevent 
any possible leakage. This-holder has the advantage that a constant check 
on the possibility of leakage can be maintained simply by covering the top 
part with water, or putting the whole in a water-bath, the temperature of 
which can be accurately measured. If the holder is of brass, troubles of 
oxidation and rusting which result in the case of iron holders can be 
avoided. 


FIG 1A CC, 


RV: 
PT, 


RV, 


PT, 


Fras. 1 anp la. 


(2) For the measurement of inlet pressure, the following set of gauges 
and manometers was provided : 

(a) Rough gauge (R@) for the rough estimation of pressure. 

(6) Calibrated Alex. Wright metal gauge (WG), for pressures from 
130 to 300 cm. Hg. , 

(c) A straight mercury manometer (MM), for pressures from 
25 to 130 cm. Hg. 

(d) A manometer (JM) with an inclination of 1 in 4, for pressures 
from 0 to 25 cm. Hg. 


Outlet pressure was measured by means of a water manometer (WM), 
separated by a calcium chloride tube (CC,) from the sample holder. 

(3) Up to rates of flow of 1 c.c./sec. a calibrated flowmeter (FM) was 
used, consisting of a resistance to flow (r), represented by a capillary of 
0-1 em. diam. and an inclined (1 in 10) toluene manometer. 

Higher rates of flow were measured by means of a Parkinson gasmeter 
(PGM) and a stopwatch. The gasmeter gave one revolution for 0-5 1. gas. 
it was separated from the sample holder by a calcium chloride tube (CC3). 

To make the measurements on the glass filter the whole sample holder 
was replaced by the glass filter which was connected to the apparatus 
in an appropriate way. . 

CALIBRATIONS. 


(1) The Alex. Wright metal gauge was calibrated against a mercury 
manometer and the calibration checked for hysteresis. 
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All the readings were converted into true values by means of a chart. 
(2) The Parkinson gasmeter was found to read correctly for 2 1. of air, 


which were collected under water, and for 4 1. which were displaced from a Rate 
standard flask by water. to Q, fi 
(3) Calibration of the flowmeter was carried out for the two purposes : Whe 
(a) To ascertain that the flow of air through the flowmeter within [f *PPT°F 
the required range was viscous, as already calculated by ascertaining on 
va 
FLOWMETER CALIBRATION © 
press 
WITH GASMETER 
OF © WOCCS BURETTE — 
0 
Table 
ments 
S one of 
from 
Q 9 A sample 
{ VA 
4 
O 
Q 
. cm. ig 
198-0 
188-4 
165-4 
144-4 
FLOWMETER READING ie 
oO 2 30 40 50 60 | 


Fia. 2. 


the Reynolds’ number for the maximum rate of flow to be measured by 
the flowmeter. 

(6) To obtain the equivalence of the flowmeter readings in c.c. /sec. 
Calibration was first carried out by means of the gasmeter, and then 
by means of a 200-c.c. burette, which was filling with water from a 
constant pressure head, acting across a fine jet, while the air so dis- 
placed was passed through the flowmeter. 


At the end of the investigation the calibration was verified by means of a 
50-c.c. burette, collecting the air, issuing from the flowmeter, under water. 
In all cases the time was taken by means of a stopwatch. The calibration 
is shown graphically on Fig. 2. 
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CORRECTIONS. 


Rate of flow as measured by the apparatus was in all cases converted 
to Q, for considerations already given. 

Whenever measurements of air were taken in contact with water the 
appropriate correction for water-vapour pressure at test temperature 
was applied—i.e., the volume was reduced by the partial volume of water- 
vapour. This was considered necessary, because only dry air came in 
contact with the sample. When the dry air comes in contact with water, 
it becomes saturated with water-vapour, so that considerations of partial 
pressure and volume cannot be neglected. 


RESULTs. 


In order that recalculations, or the introduction of corrections deemed 
necessary by other workers but neglected by the authors, may be possible, 
Table II gives the complete calculations sheet for one series of measure- 
ments taken on sample 1. The next tables give only the final results for 
one of the samples. The results for the rest of the samples may be seen 
from their respective figures. (The Jena glass gave Fig. 4 and rock 
samples Figs. 5 to 13.) 


Taste II. 


Sample I. 
Series 1. Barometer: 74-955 cm. Hg. 
Viscosity of air : 0-01773 cp. 
Temp. : 12-9° C. + 0-1° C. 


Inlet Flow- P, — Pa Q rd atm. 
press., | meter | P,, atm. | Py,atm. | ¢.¢.Jsec. | c.c./sec. atm. x 10. 
em. Hg. | reading A =. 
196-5 6-05 3-5718 0-9864 0-1128 0-04870 2-5854 12-60 0-8100 
193-0 5-93 3-5258 00-9864 0-1106 0-04830 25304 12-52 0-7960 
188-4 5-74 3-4658 0-9863 0-1070 0-04730 2.4795 12-25 0-7860 
1845 5°57 3-4138 0-9863 0-1038 0-04650 2-4275 12-04 0-7610 
1815 544 3-3738 0-9863 0-1012 0-04580 2 11:87 0-7490 
165-4 4-77 3-1658 0-9862 0-0881 0-04190 2-1796 10-85 0-6820 
160-0 4°55 3-0918 0-9862 0-0850 0-04100 2-1056 10°61 0-6600 
153-8 43 3-0080 0-9862 0-0810 0-03990 2-0226 10-33 0-6340 
149-5 415 29558 0-9862 0-077 0-03870 19696 10-03 0-6170 
1444 3-98 2:8868 0-9862 0-0741 0-03780 1: 9-78 0-5960 
140-4 3-84 2-8358 0-9861 0-0717 0-03700 18497 9 0-5800 
133-4 3:58 2-7438 0-9861 0-067 0-03540 1:7577 9-16 0-5510 
1305 3-46 2-:7048 0-9861 0-0643 0-03440 17187 8-90 0-5380 
116-3 2-96 2-5178 0-9861 0-0550 0-03100 15317 8-03 0-4800 
108-8 2-74 2-4178 0-9860 0-0510 0-02960 1-4318 7-66 0-4490 
103-2 2-54 22-3448 0-9860 0-0472 0-02790 7-22 0-4250 
97-9 2-42 2-2748 0-9860 0-0450 0-02720 1-2888 7-05 0-4040 
92-7 2-24 22058 0-9860 0-0418 0-0 1-2198 6-71 0-3820 
86-0 2-09 2-1208 0-9860 0-0391 0-02480 1:1348 6-42 0-3560 
82-4 1-95 2-0688 0-9860 0-0362 0-02340 1-0828 6-06 0-3400 
75-7 1:77 1-9818 0-9859 0-0331 0-02! 0-9959 5-70 0-3120 
70-0 1-64 1-9068 0-9859 0-0307 0-02090 0-9209 5-41 0- 2880 
64-8 1-47 1-8388 0-9859 0-0275 0-01920 0-8529 4-97 0-2670 
64-5 1-46 1-8348 0-9859 0-0273 0-01910 0-8489 4°95 0- 2660 
58-1 1-32 1-7508 0-9859 0-0249 0-01790 0-7649 4-64 0- 2400 
52:8 1-19 1-6808 0-9859 0-0222 0-01640 0-6949 4°26 0-2180 
47-6 105 1-6123 0-9859 0-0197 0-01500 0-6264 3-88 0-1965 
42-4 oo 1-5443 0-9859 0-0178 0-01380 0-5584 3-58 0-17 
37-2 0-84 1-4733 0-9858 0-0158 0-01270 |* 0-4895 3-29 0-15: 
33-4 0-74 1-4253 0-9858 0-0141 0-01160 0-4395 3-01 0-1380 
28-8 0-64 1-3648 0-9858 0-0120 0-01000 0-3790 2-59 0-1190 
25-2 0-59 1-3173 0-9858 00111 0- 0-3315 2-46 0-1040 
20-9 0-44 1-2608 * 0-9858 0-0083 0-00730 0-2750 1-89 0-0862 
14-9 0-29 1-1819 0-9858 0-0058 0-00528 0-1961 1:37 0-0615 
12-7 0-26 1-1530 0-9858 0-0050 0-00461 0-1672 1-19 0-0524 
10-7 0-23 1-:1267 0-9858 0-0044 0-00411 0-1409 1-06 0-0441 


199 
of alr, 
rom a 
es : 
rithin 
ining 
ie ‘i 
by 
en 
1a 
is- 
or. 
on 


GRUNBERG AND NISSAN: THE PERMEABILITY OF 


THe Fiow or Aqueous So.vrtions. 
The Jena glass filter was used for two purposes : 


(1) To study the flow of aqueous solutions and to establish if possible 
the relationship between the flow of gases and that of liquids. 

(2) To measure the average pore diameter and the number of pores, 
so as to assign the correct value for the parameter of length in the 
criteria governing flow. 


The apparatus used is shown on Fig. 3. It represents an ordinary flow. 
measuring apparatus, and is characterized bya constant-temperature control, 
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Fie. 3. 


Four solutions were tested: (1) distilled water, (2) 2 per cent. n-amy! 
alcohol solution, (3) two sodium chloride solutions (0-960 and 0-614 N). 
The choice of these solutions was made so as to observe the influence of 
surface forces on the flow of liquids through porous media. Fig. 14 shows 
that the flow under test was non-turbulent. A certain deviation of in- 
dividual results from the straight-line law can be observed. 


Measurement of Pore Diameter and Calculation of the Number of Pores of the 
Jena Glass Filter. , 


The measurement of the pore radii was made according to the “ bubble- 
pressure method’ as developed by Cantor, Bechold, and others. This 
method is based on the following principle : 

If a liquid, A, filling a capillary is subjected to the pressure of air, of a gas 
or a liquid B, immiscible with A, then the pressure required to displace 
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liquid A from the capillary must be equal or slightly higher than the pressure 
caused by the interfacial tension of the two phasesacross the curved meniscus. 
There is a definite relationship between the pressure required, p, the radius 
of the capillary, R, and the interfacial tension, +. This relationship may 
be expressed by : 


26 
Pp 

The Jena filter was saturated with a liquid of known surface tension, 
and pressure by air was then applied. The break-through pressure was 


7 


R= (10) 


don 
ALISODSIA NV3N 
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O12 ©: 
PRESSURE DROP ATM. 


Fie. 4. 
FLOW OF AIR THROUGH JENA GLASS FILTER. 


observed, together with the number of capillaries through which flow of air 
took place—i.e., capillaries having approximately equal pore radii. The 
experiments revealed that the variation in pore diameter was very small. 
All except thirty-two pores had radii between 1-0 to 1-2 x 10% cm. 
Substituting a value of 10-* for the pore radius in Poisseiulle’s formula 
(equation 5), n, the number of capillaries, could be ascertained, and it was 
found to be 1-4 x 105 as an average of several experiments using different 
liquids. (This gives the equivalent circular pores. The assumption of 
circularity constitutes a severe weakness in the calculations but is made 
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for expediency, as (1) no other method is known, and (2) the calculations 
are of secondary rather than primary importance in the present work.) 
The quantity flowing through a single capillary, and hence the actual 
(not the macroscopic) mean linear velocity in the different experimental 
results for the flow of air through the Jena glass filter, could be calculated. 


1 t T 
SERIES | SERIES 2 
BAROMETER.CM.Hg 74955 74-960 


TEMPERATURE “C [129201 [140200 
VISCOSITY OF AIR | O'0I7734| OOI779¢4 
SAIVPLE | 


A P 
Oo 03 O04 O05 O6 O8 
Fie. 5. 


DIscUSSION ON THE FLOow oF AIR. 


In all cases the plot of rate of flow vs. pressure drop yielded curves passing 
through the origin, but not in a single case could the test results be repre- 
sented by a single straight line passing through the origin. As the regime 
was certainly sub-turbulent, these anomalies (on the basis of Darcy's 
law) assume significance which has to be elucidated. 

Before considering the anomalies, it is well to study possible sources of 
systematic or accidental errors which may have given rise to the apparent 
anomalies. The barometer was recently checked, and found accurate, 
and the thermometer was calibrated against N.P.L. thermometers, and found 
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accurate to +0-05° C. The viscosity of air was taken from the literature 
whilst the variation in it due to variation of temperature is seen from the 
tables to be negligible. The manometric section of the apparatus is of an 
absolute type, except the steel gauge, which was calibrated against a mercury 
manometer. That the anomalies are not due to errors in pressure measure. 
ments is shown by the fact that although the pressure range (see Table II) 
covered all three instruments, no break is apparent in the curve at the 
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point of change from one instrument to another. The discontinuity on the 
curve appears to be in the range covered by the vertical mercury mano- 
meter. The usual precautions regarding verticality of scale, cleanliness 
of mercury, constancy of zero point, etc., were, of course, taken. The last 
piece or apparatus is the flowmeter, and this was calibrated by three 
different methods, in two of which dry afr was first passed through the 
flowmeter, and later measured, either by the gasmeter or a burette, and in the 
other of which air was first measured in a burette and ther displaced through 
the meter. As the curve on Fig. 2 shows, the three calibrations agree 
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well, although the first two were made before the investigation, and the last 
after the termination of all these studies. The stopwatch was checked 
against radio signals over a period of 3 hours. Thus it is concluded that, 
whilst errors due to instrumental inexactitudes are there, these cannot be 
the cause of the apparent anomalies. 
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Another possible cause for the anomalies may be the method of calcula- 
tion and corrections introduced or neglected. Thus, the assumption of 
isothermal expansion may not be correct. It is felt that the facts that the 
temperatures of the air inside the tank and in the laboratory were sub- 
stantially the same and the small pressure drops at which the anomalies 
appear, justify such an assumption. Further, adiabatic expansion, if it 
existed, would have resulted in an apparent increase in permeability at 
higher rates of flow, instead of the decrease shown.5 Again, deviations 
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from Boyle’s law at the pressures used are negligible for air. Again, 
the broad folicies involved in the method may be questioned. Foremost 
in this section is the possible condensation of water-vapour on expansion, 
Whilst the air used was not tested for complete dryness, it was felt that 
passing it through the calcium chloride dryer, 6 cm. in diameter and 16 em. 
in length, at the slow rates involved in the regions exhibiting the anomalies, 
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SAMPLE 6 4 


Ww 
(S95) 01K%G 


AP ATM/CM. 


10. 


would completely eliminate this source of error. No access of water- 
vapour to the sample through back vapour pressure was allowed. Again, 
the possibility of plugging by solid particles was investigated by the 
several repetitions of the measurements going over the high- as well as the 
low-pressure differential range (see Table III). Constant check against the 
possibility of leakage was maintained as a matter of routine. Thus it was 
concluded that the anomalies exhibited in these tests were not fortuitious, 
but represented fundamental characteristics of the systems studied. To 
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Taste III. 
Sample 7. 

Qs — x 108) ap Q, |p. — on, 
Series 1. Barometer: 74-630 cm. Hg. 

Viscosity of air : 0-01784 ep. 
Temp. : 15-0° C. + 0-2° C. 
0-0817 0-07220 { 0-2601 “380 | 0-08590 | 0-0256 0-02460 [| 0-0828 6-940 | 0-02735 
0-0631 0-05710 | 02090 16-100 | 0-06900 | 0-0231 0-02220 | 0-0723 5-980 | 0-02385 
0-0504 0-04660 | 0-1630 | 13-130 | 0-05390 | 0-0194 0 01880 0-0617 5-300 | 0-02085 
0-0343 0-03260 | 0-1103 9-200 | 0-03650 | 0-0149 0-01450 | 0-0461 4-090 | 0-01522 
0-0237 0-02280 | 0-0750 6-430 | 0-02475 | 0-0115 0-01130 | 0-0342 3-190 | 0-01130 
0-0160 0-01560 | 0-0487 4-400 | 0-01610 | 0-00907 0-00906 | 0-0290 2 0-00958 
0-0117 0-01150 | 0-0355 3-245 | 0-01172 | 0-C0741 0-00740 | 0-0210 2-090 | 0-00604 
0-00704 | 0-00703 | 0-0184 1-982 | 0-00606 | 0-00520 | 0-00519 | 0-0157 1-460 | 0-00518 
0-00464 | 0-00463 | 0-0119 1-304 | 0-00393 | 0-00409 | 0-00409 | 0-0105 1-150 | 0-00347 
0-0683 0-06140 | 0-2207 | 17-300 | 0-07290 | 0-000925; 0-000925/; 0-0013 0-00043 
0-0640 | 0-05890 | 0-1735 | 16-600 | 0-05730 | 0-0601 0-05470 | 0-1947 | 15-430 | 0-06430 
0-0479 0-04430 | 0-1551 12:500 | 0-05120 | 0-0557 0-05100 | 0-1800 14-380 | 0-05040 
0-0468 0-04350 | 0-1498 12-260 | 0-04950 | 0-0484 0-04470 | 0-1578 12-600 | 0-05210 
0-0445 | 0-04150 | 0-1446 | 11-7 0-04770 | 0-0430 0-04010 | 0-1418 | 11-300 | 0-04686 
0-037: 0-03560 | 0-1236 10-030 | 0-04080 | 0-0328 0-03100 | 0-1134 8-750 | 0-03740 
0-0344 0-032¢9 | 0-1117 9-190 | 0-03680 | 0-0293 0-02800 | 0-0946 7-900 | 0-03120 
0-0302 0-02880 | 0-0973 84120 | 0-03210 
Series 2. Barometer: 74-630 cm. Hg. 
Viscosity of air: 0-01784 cp. 
Temp. : 15-0° C’ + 0-2° C, 
0-8510 0-5050 1:3585 | 142-40 } 0-4480 0-2050 0-1595 0-5619 45-00 , 01958 
0-8030 0-4820 1-3127 | 136-00 0-4340 0-1620 0-1310 0-4701 36-90 0-1551 
0-7460 0-4550 1-2590 | 128-30 0-4160 0-1290 0-1080 0- 30-40 0-1273 
0-6750 0-4210 11944 | 118-80 0-3950 0-0962 0-0835 0-3024 23-60 0- 1000 
0-6190 0-4060 1-1367 | 114-40 0-3750 0-7560 0-0835 1-2730 | 129-40 | 0-4210 
0-5590 0-3700 1-0770 | 104-30 0-3560 0-6800 0-4220 1-1994 | 119-00 0-3960 
0-4910 0-3260 1-0033 92-00 0-3420 0-5640 0-3640 1-0820 | 102-70 0: 3580 
0-4650 0-3120 0-9684 838-00 0-3200 0-4600 0-3090 0-9626 87-10 0-3180 
0-4040 0-2780 0-8918 79-40 0-2945 0-3660 0-2570 0-8390 72-50 0-277 
0-3820 “2 0-8579 75-00 | 0-2830 0-2830 0-2080 | 0-7054 58-60 | 0-2 
0-3260 0-2340 0-7732 66-00 02550 0- 2280 0-1740 0-6037 49-10 0- 1998 
0-2810 0-2070 0-7025 58-40 0-2320 0-1660 0-1340 0-4771 37-80 0-1574 
0-2450 | 0-6346 52:10 ' 0-2090 0-1100 0-0952 0 26-90 | 0-1020 
Series 3. Barometer: 74-375 cm. Hg. 
Viscosity of air : 0-01784 cp. 
Temp. : 15-0° C. + 0-4° C, 
2-135 0-9550 2-3392 | 269-00 0-7700 2-512 1-050 2-6191 “00 [ 08630 
2-172 0-9650 2-3592 | 272-00 0-77380 1- 0-916 2-2243 | 255-50 0-7350 
2-270 0-9900 2.4391 | 279-00 0-8050 1-681 0-82 1-9673 | 231-00 0-6495 
2-380 1-020 2-5191 | 288-00 0-8310 | 
3000 T T 
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show that porous materials of the same permeability range are yielding 
approximately the same curve, a plot was made of samples 1-5 (Fig. 15), 


lO 


Fie. 15. 


To obtain a scale which would superimpose all these curves the following 
was considered : 


Considering equal pressure drops 


then 


Sample 2 was taken as a standard, and for each point “ was reduced to 


Qu x __ Ky 
A  K sample 
The values of K were taken as the average of those corresponding to 


AP — 0:3 om. and AP = 0:8 oS. It can easily be seen that for all 
cm. cm. 


these samples a change in permeability occurs at approximately 


apace. 
cm. 


Table IV gives the necessary data for this plot. 
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Corresponding data for samples 6-9 are given in Table V. 


V. 


K at AP = 0-8,| K at AP = 0-3,| Average K, K, 
darcies x 10-5. | darcies x 10-5. K sample’ 


0-867 0-796 15-1630” 
263- 308-75 0-0391 
2272-0 2248-0 0-00537 
2600-0 2565-0 0-00471 
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An attempt was made to investigate the experimental results from the 
point of view of the mean free-path theory put forward by Klinkenberg.® 
This theory proposes that the permeability of porous media to gases is a 
linear function of the reciprocal of the mean pressure, but not a function 


SAMPLE 9 
l 24 
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|! 
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of the pressure drop. In plottihg K vs. ; a straight line should be obtained, 


which, when extrapolated to zero—i.e., infinite mean pressure—should 
give the permeability of the material to liquids. This method of plotting 
was adopted, but, as Figs. 16 and 17 show, no straight line was obtained. 
The curves obtained are actually the mirror images of the curves K vs. AP, 
as shown on Figs. 18 and 19. By plotting some of Klinkenberg’s results 
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as K vs. P, — Py (i.e., the equivalent of AP (Fig. 20)), it was observed 
that there was some relationship between the variation of perme ability 
and the pressure drop. For very small values of pressure drop an increase 
of permeability with pressure drop becomes evident, while at greater values 
a straight line is obtained, which straight line may correspond to the 
assymptotic part of the curves obtained in the present work. 

The tests on the Jena glass, which was of the same type as those used by 


Klinkenberg, did not result in a definite curve when K vs. ; were plotted. 


Under certain conditions the permeability of air was lower than that to 


KLINKENBERG'S 
EXPERIMENTS ON 
FILTER. AG 
OP =O-5 ATM 
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AP =5-IOATM. _| 
+P =10- 20ATM. 
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liquids, which is difficult to explain on the slip theory connecting permeability 


and the mean pressure. 
To bring out the true relationships in these contradictory results, it was 


considered necessary to apply dimensional analysis to the various factors, 
which might be fundamental. Thus: 


S(K, P, AP, v,D,e,n)=0. . 


where D = average diameter of the capillaries ; 
e = density of the gas; 
while the notation of the other factors remains as previously given. The 
different dimensionless functions were determined and tested for validity. 
Buckingham’s x Theorem 7 had to be employed, as the number of factors 
exceeded five. The functions obtained were : 
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Reynolds’ criterion, which has always been considered as one 
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of the chief controlling factors. 


Hence the equation becomes : 


where C and C’ are constants. 
Considering the density and pressure at which volume measurements 


were taken : 
= density of air at atmospheric pressure ; 
atmospheric pressure ; 


va = velocity at atmospheric pressure ; 
= volume quantity flowing per unit time at atmospheric pressure, 


Ge). ( 
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, which is in accordance with the definition of the coefficient of 
permeability ; 


functions which will have to be tested on the experimental 
results ; 


eaP 
Pa 


paP 


Pa 


va® Pa® 


Pp? 


which can be written as : 


since 


If PY oa ) and 


constant and 
AP P 


} 


ga) 
2 — constant x va?. 


) affect the value of K fundamentally, they 


will give definite curves a plotting K vs. function, while if this is not the 
case, a spread of the points will be observed. In other words, the true rela- 
-tionship between K and P and between K and AP will be shown plotting 
K vs. the two dimensionless functions obtained, and not vs. P or AP alone. 
Figs. 21, 22, and 23 show the application of these considerations to the 
experimental results of samples 2,7,and 8. The functions are not shown for 


the other samples, as the three samples shown typify the others. 
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Sample 2 (Fig. 21) shows that the permeability coefficient K is a function 
AP x P\ 
over a limited range only. The two functions where any curve is obtainable 


26 


of a over the whole range of the curve, while X is a function of ( 
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x 
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Figs. 214 anv 21s. 


are similar. Where there is spread in the points, the accuracy of measure- 
ments was lowest. 

Sample 7 (Fig. 29) shows K as a function of as and of <a 4 over part 
of the range of observations. Again there is a portion of the plot where no 
curve can be drawn, and this portion is of lowest experimental accuracy. 
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Sample 8 (Figs. 23aand 238) shows K as a function of as and o nd 


over the whole range of observations. 

At first sight these results may appear inconsistent, because all the samples 
being porous media, should show a similar variation of the different factors, 
differing only in degree, but not in fundamentals. It is quite possible, 
however, that samples 2 and 7 are showing only sections of the genera] 
curve exemplified by sample 8. The other samples fall into line with these 
three. There are still minor differences in the general shape of the curves 
to be explained. 
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Porous media, whatever their chemical properties, have been assumed to 
behave similarly, so long as the shape of the channels through which flow is 
taking place does not alter with the flow. It is to be remembered that 
porous media are made up of different chemical substances, and should 
therefore be suspected of influencing flow in different ways. Activated 
charcoal, for instance, would profoundly affect the hydrodynamical laws 
of flow when used merely as a porous conduit. Other substances would 
show different degrees of influence. They adsorb the flowing fluid to different 
extents, and even the effective diameter in which flow would take place 
would not be the same for different conditions of pressure, temperature, ete. 

M. W. Capek ® investigated the amount of air adsorbed on Al,0,, 
Fe,0;, SiO,, and other materials which frequently enter into the com- 
position of rocks, and found that the amount of air adsorbed is quite 
considerable, 
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Adsorption of gases depends on pressure. It is therefore possible that 
adsorptive effects interfere with normal hydrodynamic flow, giving a 


variation of the coefficient of permeability with the functions ¢, (Ga) 


and 6( 7), or rather with the dimensionless and more fundamental 


sp 


criteria and (45). Finally, the third function (.2.) was con- 


ou vDe 
sidered. This is Reynolds’ criterion, which is believed to determine the 
nature of flow in porous media under different conditions. It is known 


that within the ordinary range of the viscous region K is inyariant with 


+), As the flow was very slow, however, it is pertinent to assume that 
K is a function of (+). 


K = = * 


Since = is constant, it may be written : 


Even if this 


It was impossible to evaluate D for the rock samples. 
should have been possible by the method used for the Jena glass filter, the 
variation in pore size would have been too great to allow any conclusion 
to be reached. It was therefore thought advisable to evaluate the correct 
values of 

vDe 


for the Jena glass filter plotting it vs.K. This gave a curve showing a 
maximum. It was then assumed that an average value of Dexisted forevery 
sample and that this value was constant for every sample. Thus, plotting 
Qa vs. K for every sample, curves were obtained showing a minimum and a 
maximum. It was then assumed that the maximum obtained with the rock 
samples occurred at the same value of Reynolds’ number as that for the 
glass filter. Thus the permeability function of the rock samples were 
plotted on the same chart as that of the rock samples, so that the maxima 
were below each other. Fig. 24 shows the similarity of the curves thus 
obtained. 

It may therefore be assumed that a general curve exists on which the 
coefficient of permeability drops to a inieun at a value of Reynolds’ 
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criterion of about 0-2, rising to a maximum at a value of Reynolds’ 
criterion between 1 and 2 (taking the parameter for length the mean dia meter 
of the pores). For higher values of Reynolds’ criterion the value of K 
becomes more and more independent of the criterion, assuming a constant 
value, which would correspond to true viscous flow. It may therefore be 
concluded that there is the probability of the existence of a regime of 
subviscous flow for gases, where the relationship between the pressure drop 
and the velocity of flow is not governed by the simple linear equation of 
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viscous flow. This fact is also evident from Figs. 4 to 13 directly, as 
discussed above. 

An attempt was made to obtain some kind of mathematical formulation 
connecting the pressure drop with the velocity of flow. 

An increase of permeability could only be expressed by the presence of 8 
fractional power of v in the permeability equation. Examining the lower 
portions of the curves of samples 7-9, on which the increase of permeability 
is quite conspicuous, they were found to approximate to 
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The equation for the whole curve was then determined by graphical 
methods, and found to approximate to 


. . .. (18) 


An attempt to apply this equation to samples 1-6 failed. It is, however, 
important to mention that these also failed to obey equation 9, i.e. : 


= av + br?. 

The fact that maxima and minima are obtained for small values of 
Reynolds’ number is of importance for pressure conservation in gas-wells, 
as the ratio of yield to pressure drop may be as much as 100 per cent. 
higher at the maximum permeability than for slightly lower or appreciably 
higher values of Reynolds’ number. Calculations show that the linear 
velocity in the formation around gas-wells is very small. In a gas-well 
producing one million cubic feet per day measured at atmospheric pressure, 
from a stratum 30 feet thick with a borehole of 6 inches diameter; the linear 
velocity of the gas at the sand face is as low as 3-75 cm. per second at 
atmospheric pressure. Ata distance of 2 feet from the centre of the well the 
velocity is 0-94 em. per second, and at 10 feet from the centre of the well 
the velocity is 0-19 cm. per second. At the higher pressures prevailing 
in the formation the actual mean linear speeds are smaller than these values. 
Considering the small pore diameters in the formation and the small 
linear velocity, the values of Reynolds’ numbers under which production is 
taking place are very small. As the actual calculation of Reynolds’ number 
presents some difficulftes, it may be sufficient to use the volumetric rate per 
unit area of sand face and at prevailing pressures to obtain a value of 
permeability in tests designed to determine optimum rates of flow. 

To summarize: For higher values of Reynolds’ number, but below 
turbulence, the permeability of porous media to gases assumes a constant 
value which is in accordance with the conditions of viscous flow. For low 
values of Reynolds’ number the existence of a sub-viscous regime of flow . 
for gases must be assumed. It is probable that flow in gas reservoirs is 
ordinarily in the sub-viscous regime for the major portion of the reservoir. 
Thus tests should be performed with due regard to the regime obtaining 
in the reservoir. 


Discussion on the Flow of Aqueous Solutions through the Jena Glass Filter. 


The flow of aqueous solutions showed that the existence of the sub-viscoys 
regime could not be established in the case of liquids. The curve Q vs. P 
was a straight line, although some deviation could be observed for each 
solution. There remained to be found whether the adsorptive effects 
which had been assumed in the case of gases existed for liquids. 

The coefficients of permeability were different for the different liquids ; 
further, for each liquid, K showed a variation with temperature. This fact 
is brought out by Fig. 25, on which the coefficient of permeability to different 
solutions is plotted against the temperature. It can be seen from the 
graph that every solution gave a straight line with a definite slope to the 
horizontal. As the viscosity had already been considered in calculating 
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the coefficient of permeability, it was postulated that the viscosity was not 
the only property influencing flow. 

The conclusion reached was that the effective cross-section under viscous 
flow was different for different liquids, due to differences in surface energy 
and consequent differences in thickness of adsorbed layers. As the pressure 
drop varied with a high power of the radius, small differences in the thick- 
ness of these layers could account for the variation in K observed. 

A physical interpretation of this fact may be obtained by dimensional 
analysis. It may be assumed that instead of the viscosity being the sole 
criterion as a property of liquids governing flow, viscosity », density p, and 
the product of surface tension and the average pore diameter d, are all 
important. Then a dimensionless function obtained is : 


(ca) 


As d is a constant for each rock sample, and is unknown for different 
samples, it is omitted from the criterion. Thus: 


It is to be remembered that this function is only dimensionless when the 
product (od) is used to derive it. No dimensionless function can be obtained 
from », pe and o alone. 
It is found experimentally that : 
2 
. ail 
log K = log K, + 0-1 log (20) 


Thus the dimensions of K, are those of K. To render the equation 


2 
dimensionless in both variables, K and f (), 


log = 041 log (©) 


Fig. 26 gives the plot-for this equation, and it is seen that all four liquids 
yield points falling on a single curve. The function is in reality 
l 

and absolute viscosity to that due to surface-energy effects. K is then a 
measure of permeability as governed mostly by viscosity phenomena, and 
K, of permeability as governed by surface energy. In other words, where 


), and thus may be taken as the ratio of resistance due to kinematic 


pand . are small and c is large, flow is affected to agreater extent by surface- 


energy effects than in cases where » and “ are large and oc is small. This 


is in accord with ordinary experience. 

It was concluded that the coefficient of permeability assumes different 
values with different liquids, and is also dependent on the temperature for 
each liquid. The slope of the graphs of K vs. temperature in degrees 
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Centigrade is the same for all the li 
Thus : 


quids tested and is equal to 0-80 md/°¢. 


k=a—b(#C.) 
where a is a characteristic constant of the liquid used. It is to be noted 
that if this equation is the same for all liquids and media used, then from 
a determination of permeability of a medium to a particular liquid at a single 
1-200 


+ 


(toc. scace ) 


temperature the permeability at other temperatures can be calculated. 
, Further work is necessary, however, before this condition can be assumed 
to be operating. 


CONCLUSIONS AND RECOMMENDATIONS. 


(1) The permeability of porous media to gases and liquids presents two 
different problems and in the present work no correlation between these two 
is apparent. 
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(2) Permeability to gases is dependent mainly on Reynolds’ number, 
which at small values gives rise to maxima and minima. For higher 
values of Reynolds’ number, but preceding the turbulent regime, per- 
meability assumes a constant value corresponding to viscous conditions. 
Thus for low rates of gaseous flow in rocks, tests should be carried out at 
Reynolds’ numbers corresponding to those obtaining underground. The 
“ Reynolds’ number ” may be taken simply as the volumetric rate of flow 
per unit area of sand, the area being perpendicular to the main line of flow. 
In other words, the mean linear speed of the gas is taken as the “ Reynolds’ 
number.”’ It is further recommended that tests be carried out under a 
mean pressure equal to that existing underground. 

(3) Permeability to liquids is dependent not only on the viscosity, but 
also on other properties of the liquid which cause a variation in the effective 
cross-section of the pores under viscous flow. The coefficient of permeability 
assumes different values for different liquids having different surface-tension 
values. A single curve was obtained for the different liquid used by 


2 
plotting against the criterion (©). The permeability also varies with the 


temperature for the same liquid, but a correction for this can be obtained 
for each liquid : 


K=a—bFC. 


Further work is necessary before these equations can be taken as uni- 
versal. The effects of the temperature are fully accounted for in the 


2 
criterion (« ). Thus for correlation between different liquids at different 


temperatures this criterion should be used. For correlation with the same 
liquid at different temperatures the more specific temperature equation 
may be used. 


Department of Oil Engineering and Refining, 
University of Birmingham. 
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THE PHASE TRANSFORMATIONS OF NORMAL 
PARAFFINS. 


By C. G. Gray, B.Se., A.H.-W.C., A.C. 


INTRODUCTION. 


In recent years a considerable amount of information has been published 
on the crystalline habits of normal paraffins and other long-chain com. 
pounds. These studies have shown that polymorphism is common in 
such compounds, and that crystalline transitions frequently occur in the 
solid phase. The solid phase transitions have been studied in some detail 
by various experimenters, each working with a more or less restricted 
range of compounds, but no attempt appears to have been made to cor. 
relate the results observed in these restricted zones and to form a complete 
picture of the solid phase relationships for a given class of compound. 
This task is attempted for the normal paraffins in the present paper. 
The stydy of the solid phase relationships in the normal paraffins throws 
some interesting light on a number of difficulties and apparent incon. 
sistencies in the literature on the crystallography of petroleum waxes. 
The crystallography of the pure n-paraffins is in itself a subject of some 
complexity ; it is therefore not surprising that the interpretation of results 
on commercial waxes is full of pitfalls. 


Tue CRYSTALLINE MOopIFICATIONS OF NORMAL PARAFFINS. 


Pure normal paraffins are known to crystallize in at least four distinet 
modifications :— 


(1) Hexagonal System, the hydrocarbon molecules being packed side 
by side “like hexagonal pencils ” (Miiller), the long axes being perpendi- 
cular to the planes containing the ends of the molecules. This is shown 
diagrammatically in Fig. la. The molecules rotate freely about their 
long axes. 

(2) Orthorhombic System (Fig. 16). In this modification the molecules 
are packed side by side perpendicular to the planes containing their ends, 
i.e. the (001) planes, but the packing is no longer hexagonally symmetrical. 
Two unequal side spacings are shown. The molecules do not possess 
freedom of rotation, but show libration about their mean positions. 

(3) System uncertain (monoclinic or triclinic). Two unequal side packing 
spaces are shown, and the molecules are tilted at an angle of 73° to the 
planes containing their ends (Fig. Ic). 

(4) System uncertain (monoclinic or triclinic). Two unequal side spac- 
ings are shown, and the molecules are tilted at an angle of 61° 30’ to the 
planes on which their ends are situated (Fig. 1d). 


In the literature, a great deal of confusion exists over the nomenclature 
of these four modifications. In the subsequent discussion, the following 
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nomenclature will be used throughout, and it is suggested that it would 
provide a convenient standard notation :— 


aform: Vertical rotating chains (hexagonal). 

form: Vertical non-rotating chains (orthorhombic). 
y form : chains tilted 73° to (001) planes. 

sform : chains tilted 61° 30’ to (001) planes. 


The various crystalline modifications described above are not observed 
in all normal paraffins, the grystalline habit being governed by the number 


VERTICAL LATTICE (HEXAGONAL), 


+ VERTICAL LATTICE (ORTHORHOMBIC), 
te eo! Latnee turep 73° 
e 
(SysTEM UNCERTAIN). 
Fro. 1 


CRYSTAL SYSTEMS OF NORMAL PARAFFINS. 


of carbon atoms in the molecule. The ranges of n in which each modifica- 
tion has so far been observed are as follows :— 


a: Odd paraffins from C,, upwards; even from C,, upwards. 

8: Odd from C,, upwards; even from C,, upwards. 

y: Odd from C; to C, inclusive; even from C, to C,, inclusive. 

8: Even only, C,, to C3, inclusive. This modification has been 
obtained only by crystallization from solvents. It has not been 
obtained by crystallization of melts or by transition from one of the 
other modifications. 


Tue MELTING Pornts OF n-PARAFFINS. 


The melting points of the normal paraffins below C,, show odd-even 
alternation, the values lying on two smooth curves, that of the even 
members lying above that of the odd. The two curves converge at 
n= 20. Above n = 20 the values for both odd and even members lie 
on one smooth curve. This curve has been assumed to continue up to 
the highest known paraffins, although erratic values mage been found in 


the region to Cyo. 
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THEeory OF MELTING. 


W. E. Garner and his co-workers have shown that the heats of crystalliza. 
tion of the vertical rotating («) forms of normal paraffins increase linearly 
with the number of methylene groups in the molecule. The two termina] 
methyl groups, on the other hand, contribute a negative quantity to the 
molecular heat of crystallization. The contributions to the entropy of 
crystallization are in a similar sense, and it is concluded that “ there is q 
marked tendency for these ends to enter the liquid state.”’ 

These facts throw light on the relationship between the melting point 
and the number of carbon atoms in the molecule for those normal paraffins 
of which the vertical rotating forms are stable at the melting point. 

The latent heat and entropy contributions make it clear that the bonds 
holding a paraffin molecule in the (vertical) crystal lattice are associated 
with the (n—2) methylene groups, the terminal methyl groups being 
repelled by the lattice. When the crystal melts, these (n—2) bonds have 
to be broken, the energy required being provided by (a) the potential 
energy of the terminal groups and (5) the vibrational energy of the molecule. 

According to the quantum theory of molecular heats, the vibrational 
energy of the molecule is composed of two separate contributions (a) the 
energy of vibration of the molecule as a unit and (5) the intramolecular 
energy of vibration of the atoms within the framework of the molecule. 
The molecular vibration energy is equal to RT’ per mole for each degree 
of vibrational freedom. Taking into account only the energy in the 
direction parallel to the (n—2) lattice bonds, and distributing this quantity 
RT equally over these bonds, the contribution to each bond is RT (n-2) 
per mole. 

The atomic intramolecular vibration energy is given by Einstein’s fune- 
tion hy /e 
this energy approximates closely to R(7'-C) per mole for each degree of 
freedom, where C is a constant for a given molecular structure. The total 
vibrational energy in the direction of the lattice bonds is therefore 
R(T-C) + RT/(n-2) per bond per mole. 

The fact that the heats of crystallization of the methylene groups and 
of the terminal groups respectively are constant above C,, irrespective of 
the temperature of crystallization, shows that the bond energy and the 
potential energy of the terminal groups are constant for similar crystalline 
forms, irrespective of the temperature of crystallization. 

Let the bond energy be Z, per mole and the potential energy of the 
methyl groups distributed over the (n-2) bonds, be 2£,/(n—2) per 
bond per mole. It is reasonable to assume that at melting, the bond 
energy is neutralized by the combined potential and vibrational energies, 
so that 


1, and it is known that except in the region of absolute zero 


2B, RP 


= + R(T-C) + ——. 


(n-2) (n-2) 
This equation reduces to the form . 
T = A — B/(n-1) where A, B, are constants. 


Over the range Cy, to C3,, the vertical rotating forms of the n-paraffins 
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Taste I. 

Melting Points of Normal Paraffin Hydrocarbons. 
| | | At? C. 
(cale. 
minus 
obs.). 


2143-41 6 | 

(8 form) 9 

908-86 

1908: il 

Odd: t= E1779 12 

(a form ?) 14 

15 
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78-9 
2- 
84-7 
197-4 2266-7 44 85-9 86- 
37-4 — 50 92-1 92-1, 93-0 * . 
52 93-8 94 0-2 
(8 form 54 95 0-4 
| 60 99-6 98-9 0-7 
62 | 100-8 100-5 0-3 
| 64 102-0 102 0-0 
66 | 103-1 103-6 —0-5 
67 103-6 104-1 
70 | 105-02 105-25 —0-23 
* Egloff quotes two values without discrimination. 


+ According to Egloff, “ impure.” 
are stable at the melting point, and the melting points for this range, 
calculated from the equation T(°K) = 396-4 — i, are shown in Table I, 
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together with the observed melting points published by Egloff? The 
constants A and B were evaluated by substituting the observed values of 
T for and 

Above C3,, and below Cy9, the values of 7’ calculated from this equation 


show marked deviations from observed values. It is known that the even 


paraffins below Cy) crystallize with tilted chains, the melting point curves 
of the vertical and tilted forms intersecting at about Cy.2 The odd 
paraffins above C, do not crystallize with tilted chains; nevertheless, the 
melting point curve of the odd paraffins also shows a discontinuity at 
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Grapu 1. 
INTERSECTION OF a- AND §-MELTING POINT CURVES AT C5. 


about C5. The equations of the “ odd ” and “ even ” curves in this zone 
have been calculated from the observed values, and it has been found 
that these curves are give by the equations :— 


1908-86 

n +- 1-779 
2143-41 

n + 3-7069 


T(°K) = 396-77 — odd series, C19. 


T(°K) = 399-81 — even series Cg 99. 
The melting points calculated from these equations are compared with 


observed values in Table I. 
The data available on the paraffins with more than thirty-six carbon 
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atoms in the chain are neither so complete nor so reliable as those for the 
lower members. Nevertheless it is here suggested that the discontinuity 
at Cy, is due to a definite change in the form of crystal which is stable at 
the melting point above C,,. The evidence for this suggestion will appear 
in the next section. The melting points in the zone above C,, lie on the 


curve: (°C. = 137-4 — — as is shown in Table I and in Graph 1. 


Tue PHase TRANSITIONS. 


The curve for melting point above C,, depressed through 2-4° C., con- 
tinues as a transition in the solid phase for paraffins below C,,. The 
“ first transition points on cooling ” reported by Piper and his associates ¢ 
fot Cyg—Cy9, Cgy-Cyg are plotted in Graph 1, and are compared in 
Table II with values calculated from the equation 


£0. = 135 — 


Piper and his co-workers confirmed that for paraffins above C,, this 
transition was from the vertical rotating form to the vertical non-rotating 
form. The tilted forms (angle of tilt 61° 30’) obtained by crystallizing 
the even paraffins in this range from solvents pass monotropically to the 
vertical rotating forms, the temperatures at which this takes place laying 
0-5-1-5° C. above the “ First transition on cooling” line, except for Cy, 
and C,,, which show a difference of about 3-4° C. between the two transitions. 


Taste II. 
Transition Temperatures. 


Piper's Calculated transition At 
“ transition point (a —> 8) observed 
point, = 135-0 — 22067 | minus 
cooling.”’ n calculated. 
26 48-3 47-8 
| 27 51-0 51-0 0-0 
| 28 54-0 54-1 
29 55-8 56-8 —1-0 
| 30 58-0 59-4 —1-4 
| 61-8 61-9 
| 32 63-9 64-2 —0-3 
| 4 68-5 68-3 +0-2 
| 35 70-5 70-2 +03 
36 72:5 72-0 


Piper’s “ first transition point on cooling,” with the nomenclature 
adopted in the section on Crystalline Modifications, is the «<—>§ transition 
for both odd and even series. The continuity of this «<—+>8 transition line 
with®the melting-point line for paraffins above C,, is a strong indication 
that the $ form crystallizes from the melt above C,,, and not the « form, 
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as in the Cy), range. Piper’s “transition on heating” at which the 
8 form passes to « lies a little above the 8 setting point line in the solid 
phase, but King and Garner ' have pointed out that for the longer chaing 
the vertical forms are probably stable at the melting point. 

The identification of Piper’s “cooling” and “ heating ”’ transitions 
with the a<—+8 (odd-ev€n) and 3——> «a (even) changes, respectively, 
immediately raises the question of the nature of the transition line which 
runs roughly parallel to and some 5° C. below the « melting point curve 
in the range C,,~C,., as this line has been assumed to be the line of transi. 
tion from « to 8 in this range.6 This transition has apparently been 
noted in even paraffins only, no data on odd paraffins in the range having 
been traced. 

Kolvoort > and Miiller® examined the pliase transitions of n-C,,H,,, 
and found that the two transitions, which take place approximately 5° 
and 10° C. respectively below the melting point, are reversible. On the 
other hand, the upper transition point in C,, and C,, is monotropic, as 
shown by Piper. Above C,, this transition does not appear at all. Piper 
and his associates regarded their sample of triacontane as anomalous in 
this respect, and continued the curve of the upper transition for C,, and 
C,, to join up with those of C,, and C,,. The evidence available, however, 
indicates that the upper transition below C,, does not correspond with 
the upper transition in C5, and higher hydrocarbons. This is shown in 
Fig. 2, which shows the observed transition points over the range C,,-C,,, 
and the suggested scheme of correspondence between them. 

It is perhaps significant that the melting-point curve for C,, and over 
cuts the line of upper transition temperatures for C,, and under at the point 
corresponding to C9, at which point the upper transition (heating) is 
suddenly depressed (Piper), and is thereafter identified as the (8 —-> 2) 
transition. 

The y melting-point curve for even paraffins intersects the « melting- 
point curve at C,, approximately,® and it is suggested here that the upper 
transition in the range C,)—C,, is the continuation of this line as a transition 
a<—>y for even members. Data on odd members do not appear to be 
available. The odd members above Cy, however, are not known to exist 
in a tilted form. 

This interpretation is at variance with the conclusions of Kolvoort,® 
who classified the upper transition for n—C,, as the point of (enantiotropic) 
conversion from hexagonal « to orthorhombic 8. Kolvoort did not, 
however, claim to have established the crystalline systems in which his 
preparations separated. His observations were confined to the angles 
between the lines of intersection of the (001) and (hk0) planes, and were 
thus incapable of differentiating between vertical and tilted forms. It 
was assumed that Miiller’s X-ray data on C,, would also be applicable to 
C.,. This assumption, however, is no longer justified. 

Miiller showed that the lower transition in C,, corresponded with the 
sudden displacement of the (200) reflection in the X-ray spectrum (40- 
41° C.). Thereafter, with rising temperature, the (200) spacing increased 
continuously until it became identical with the (110) spacing (46° C.), 
this being interpreted as a continuous transition of an orthorhombic 
structure into a radially symmetrical (hexagonal) form in which the chains 
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are free to rotate. It is, however, possible that the lower transition is 
from one tilted modification to another (e.g. triclinic —> monoclinic), 
and that the chains do not assume a position vertical to (001) until the 
upper transition point is reached, at which point the rotation sets in and 
the molecules take up a vertical position. This point can only be settled 
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PHASE TRANSFORMATIONS OF NORMAL PARAFFINS. 
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3: B m.pt. (odd) (B >a). 
6 m.pt. (even) (monotropic) (8 -> a). 
4: B set pt. (odd-even) (a > 
5: y m.pt. (enantrotropic below C,,) (even) (a <—> y). 
6: 5 m.pt. (enantrotropic below C,,) (even) (y<—>8). 
7: y m.pt. (even). 
8: am.pt. (odd). 


by measurement of the long spacings over the temperature range involved, 
evidence from short spacings being insufficient to settle the question 
definitely. 

Yannaquis 7 and Bijvoet® concluded from X-ray data (the latter on 
pure n-C,,) that the upper transition point in this range corresponds with 
a change from a hexagonal « form to a tilted form. According to Yan- 
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naquis the tilted form is triclinic; according to Bijvoet it is monoclinic 
for C,,. Bijvoet gives the following scheme for C,,. 


Hexagonal 

46-5° C. 
Monoclinic 1 

42° C. 


Monoclinic 2 (metastable) 
or Triclinic (stable at room temperatures) 


This agrees with the hypothesis that the upper transition for even paraffins 
in this zone is the «<—»y transition. 


SUMMARY AND CONCLUSIONS. 


The equilibrium diagram shown in Fig. 2 is advanced as being in ac. 
cordance with the available data on the melting and solid-phase transitions 
of normal paraffins. 

Theoretical considerations lead to the formula 7 = A — B/{n—1) for the 
relationship between melting point and number of carbon atoms in those 
n-paraffins in which the « form is stable at the melting point. The formula 
is applicable over the range n = 20-36 inclusive, and is in excellent agree. 
ment with observed values. 

Above C,, the relationship between melting point and number of carbon 
atoms is of a different type, and the curve for melting point in this zone 
continues as a transition («<—>8) in the solid phase below C;,, suggesting 
that the 8 form is stable at the melting point above C,,. 

The upper transition line in the zone C,.—C,, is probably the continuation, 
as a transition («<—»y) in the solid phase, of the y melting-point curve for 
even paraffins below Cy». 
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CALIBRATION OF C.F.R. REFERENCE FUELS. 


Calibration of Secondary Reference Fuel F.4+- 4 mls. T.E.L./Imp. Gln. 
and C12 +4 mls. T.E.L./Imp. Gln. by C.F.R. Motor Method I.P. 
144/42(7’) and the 17° Motor Method I.P. 43/42(T)). 

Report of Sub-Committee No. 5—Engine Tests. 
Aviation Fuel Knock-Rating Panel. 


INTRODUCTION. 


Ix accordance with standing procedure, the above calibration has been 
considered on the introduction of F.4 Reference Fuel. As F.4 appeared 
to be similar to F.3 it was-decided to determine whether the calibration 
for leaded F.3—C.12 would be accurate enough for use with the latest batch. 


TEsTs. 
Blends consisting of (50% F.4 + 50% C.12) + 4 mls. T.E.L./1.G. 
and (90% F.4 + 10% C.12) + 4 mls. T.E.L./LG. 


were rated in terms of leaded F.3 and C.12 in ten C.F.R. engines by 17° 
Motor Method. The results of test are given in Table I. 


I. 


(50% F.4 in C.12) + 4 mls. (90% F.4 in C.12) + 4 mls. 
C.F.R. Engine. | T.E.L./I.G. in (% F.3 in C._?) | T.E.L./1L.G. in (% F.3 in C.12) 
+ 4 mls. T.E.L./1.G. + 4 mls. T.E.L./L.G. 
Diff. Diff. ' 
A 49-8 —0-2 89-6 —0-6 
B 91-0 +0°8 
Cc 52-2 * +2-2 93-0 * +2°8 
D 50-0 0 90-0 —0-2 
E 50-4. +0-4 90-0 —0-2 
F 50-3 90-3 +01 
G 49-7 —0-3 90-0 —0-2 
H 49-8 —0-2 90-2 0 
I 50-0 0 90-0 —0-2 ! 
J 50-0 0 91-0 +0°8 
Average 50-0 90-2 


* Ratings on Engine ‘‘C”’ have not been included in the average, as they were 
obviously in error. 


CONCLUSION. 


As the average difference is negligible (less than 0-1 O.N.), it is recom- 
mended that the calibration tables issued by the Institute of Petroleum 
for F.3+4 mls. T.E.L./L.G. in C.12+ 4 mls. T.E.L./L.G.* shall be 
employed for leaded F.4 and C.12. 


* J. Inst. Petrol., 1941, 27 (211), 192. 
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The amended Table is reproduced herewith as Table ITI. 


II, 
Octane-Number Calibration Tables C.F.R. Motor Method (I.P. 44/42(T)) up to 1 
O.N., thence by “17° Motor Method” (I.P. 43/42(7')). 
F.4 + 4 mls, T.E.L./LG. in C.12 + 4 mls. T.E.L./1L.G. 


% F.4 + T.E.L. N. % F.4 + T.E.L. OLN. 


103-0 
103-4 
103-9 
104-4 
104-8 


© 


ome 


105-2 
105-7 
106-2 
106-6 
107-0 


107-5 
108-0 
108-5 
109-0 
109-5 


110-0 
110-5 
111-0 
111-5 
112-1 


112-7 
113-3 
114-0 
114-6 
115-3 


116-0 


0 50 
4 54 
6 56 
8 58 
10 60 
12 62 
14 64 
16 66 
18 68 
20 70 
22 72 
26 76 ig 
238 78 
30 80 
32 82 
34 
36 l 86 + 
38 1 88 ree 
40 101-0 90 ae 
42 101-4 92 at 
44 101-8 94 
46 102-2 96 
102-6 98 
100 


= 


